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Abstract  A  cold  thermoregulatory  model  (CTM)  was 
applied  to  data  from  partially  immersed  subjects 
divided  into  normal  (NF)  or  low  fat  (LF)  groups  in 
order  to  validate  CTM  during  immersion  at  two  depths 
and  to  examine  mechanisms  underlying  the  individual 
differences.  CTM  defines  thermal  characteristics,  e.g. 
surface  area  and  maximal  shivering  intensity,  using 
height,  weight,  fat  %,  age  and  V02^^.  Ten  clothed 
subjects,  5  NF  (15-19%)  and  5  LF  (8.1-14.7%),  were 
immersed  in  both  10  and  15°C  water  at  chest  (CH)  and 
waist  (WA)  level.  Environmental  and  clothing  inputs 
for  CTM  were  weighted  to  adjust  for  the  ratio  of  skin 
surface  area  covered  by  either  air  or  water  at  various 
immersion  depths.  Predicted  core  temperature  (Tj.) 
responses  for  each  individual  trial  were  compared  with 
measured  data.  There  were  no  significant  differences 
(F  >  0.05)  between  measured  and  predicted  for 
NF  at  aU  four  conditions.  In  contrast,  for  the  LF  group, 
the  predicted  responses  were  all  higher  than  mea¬ 
sured  (P  <  0.05).  However,  predicted  agreed  closer 
with  measured  for  LF  when  leg  muscle  blood  flow 
was  increased  in  the  simulation.  This  suggests  that 
blood  flow  may  contribute  to  the  rapid  decline  in 
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observed  in  LF  and  its  variance  may  cause  in  part  the 
individual  differences  in  responses.  CTM  predicts 
responses  to  immersion  at  various  depths  with  accept¬ 
able  accuracy  for  NF  individuals  in  this  study  and  can 
be  adapted  to  non-uniform  environments. 

Keywords  Model  •  Cold  stress  •  Hypothermia  • 
Thermoregulation  ■  Muscle  blood  flow 


Introduction 

One  of  the  difficulties  in  establishing  criteria  for  cold 
survival  or  injury  prevention  and  in  developing  cold 
thermoregulatory  models  is  the  intra-and  inter-individ¬ 
ual  differences  in  the  cold  response.  For  example,  in 
Molnar’s  survival  time  curve,  survival  times  ranged 
from  2  to  50  h  during  immersion  in  ocean  water  of  20- 
24°C  (Molnar  1946).  There  is  no  simple  explanation  as 
to  why  such  large  inter-individual  differences  exist  and 
what  the  parameters  are  in  determining  individual 
responses.  Potential  factors  include  height,  weight, 
body  fat,  shivering  capacity,  blood  flow  and  immersion 
depth,  etc.  These  factors  impact  either  physiological 
(shivering,  vasoconstriction,  vasodilation)  or  biophysi¬ 
cal  (heat  transfer  from  the  core  to  the  body  surface 
through  conduction  and  blood  flow  convection)  pro¬ 
cesses  of  heat  balance. 

Body  fat  is  a  significant  factor  in  the  relationship 
between  human  thermal  equilibrium  and  environmen¬ 
tal  cold  stress  (Toner  and  McArdle  1996;  Stocks  et  al. 
2004).  The  fat  influences  the  physiological  and  internal 
thermal  processes  in  two  ways:  (1)  it  provides  insula¬ 
tion  since  its  conductivity  is  low  (Werner  and  Buse 
1988;  Toner  and  McArdle  1988);  and  (2)  it  is  typically 
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associated  with  low  metabolism  and  thus  attenuates 
shivering  heat  production  (Tikuisis  et  al.  1988).  Most 
studies  report  a  inverse  relation  between  body  fat  com¬ 
position  and  cooling  rate  (Hayward  and  Keatinge 
1981;  Veicsteinas  etal.  1982;  McArdle  etal.  1984) 
although  one  study  showed  that  subjects  with  low  fat 
maintained  a  similar  core  temperature  (7^).  compared 
to  their  high  fat  counterparts,  due  to  a  significantly 
greater  shivering  thermogenesis  during  immersion  in 
18,  22  or  26°C  water  (Glickman-Weiss  etal.  1991). 
Since  body  fat  plays  an  important  role  during  cold 
exposure,  it  is  necessary  for  thermoregulatory  models 
to  take  this  factor  into  account  in  order  to  differentiate 

responses  in  populations  with  various  body  fat  com¬ 
positions.  Within  population,  body  fat  composition  can 
be  changed  significantly  after  intensive  training  (Friedl 
et  al.  1994)  and  consequently  the  capability  in  defend¬ 
ing  cold  may  change  too. 

Muscle  blood  flow  is  another  factor  which  influences 
heat  balance  and  internal  heat  transfer  process.  How¬ 
ever,  quantitative  information  on  relationships 
between  muscle  blood  flow  and  human  thermal 
responses  to  cold  is  limited,  because  muscle  blood  flow 
is  difficult  to  measure  accurately.  A  recent  study  mea¬ 
sured  muscle  blood  flow  using  an  optic  fiber  and  dem¬ 
onstrated  the  potential  contribution  of  muscle  blood 
flow  to  changes  in  after  exercise  in  25°C  air  (Kenny 
etal.  2006).  The  contribution  of  muscle  blood  flow 
would  be  expected  to  be  greater  during  cold  water 
immersion  due  to  higher  temperature  gradients 
between  the  core  and  extremities.  One  optional 
approach  to  study  the  impact  of  muscle  blood  flow  on 
heat  balance  and  responses  could  be  through  simu¬ 
lation  using  human  thermoregulatory  models. 

Biophysical  analysis  of  heat  balance  becomes  more 
complicated  during  partial  immersion  since  the  exter¬ 
nal  environment  is  not  homogenous,  i.e.  people  are 
simultaneously  exposed  to  both  air  and  water  environ¬ 
ments.  Consequently,  no  cold  model  has  validated  its 
prediction  capacity  under  partial  immersion  condi¬ 
tions,  A  study  was  conducted  to  examine  the  role  of 
immersion  level  and  water  temperature  on  during 
light  exercise.  The  subjects  were  ten  young  healthy 
men  with  body  fat  compositions  ranging  from  8.1  to 
19%.  The  present  study  applied  our  cold  thermoregu¬ 
latory  model  (CTM)  (Xu  et  al.  2005)  to  analyze  this 
physiological  data.  CTM  uses  five  biophysical  charac¬ 
teristics,  including  the  percentage  of  body  fat,  to  define 
the  overall  thermal  characteristics/properties  of  each 
individual.  The  goals  were  (1)  to  determine  whether 
CTM  can  predict  thermal  responses  for  each  individual 
with  various  fat  composition  during  immersion  in  cold 
water  at  two  depths;  and  (2)  to  find  out  which  factors 


or  mechanisms  (i.e.  fat  %,  shivering  or  muscle  blood 
flow)  may  contribute  to  large  individual  differences, 
especially  in  individuals  with  low  body  fat  composition. 
By  modeling  the  effects  of  different  immersion  levels  in 
combination  with  body  composition,  the  potential  life- 
threatening  hazard  of  partial  immersion  may  be  better 
understood,  and  more  active  preventive  measures  can 
be  implemented  or  public  warnings  issued. 

Methods 

Cold  thermoregulatory  model 

Cold  thermoregulatory  model  was  derived  from  the 
previous  version  of  a  thermoregulatory  model  devel¬ 
oped  by  Werner  and  Webb  (1993),  which  was  based  on 
the  essential  developments  of  Stolwijk  and  Hardy 
(1977).  It  was  a  six-cylinder  model,  with  each  cylinder 
consisting  of  a  core  and  shell  layer.  Recently,  we  added 
muscle,  fat  and  clothing  layers  (Xu  and  Werner  1997) 
and  incorporated  a  conceptual  model  for  shivering 
intensity  and  fatigue  into  the  current  CTM  configura¬ 
tion  (Xu  et  al.  2005).  This  has  improved  its  perfor¬ 
mance  for  predicting  human  responses  to  long-term 
cold  exposure.  Each  cylinder  is  now  divided  into  con¬ 
centric  compartments  representing  the  core,  muscle,  fat, 
and  skin.  The  outer  cylinder  has  an  additional  clothing 
layer.  Circulation  is  represented  as  a  one-loop  circula¬ 
tory  system  and  is  an  independent  compartment.  Thus, 
the  human  body  is  represented  by  25  compartments. 
The  sizes  of  the  compartments  are  determined  from 
height,  weight  and  percentage  of  body  fat  (Xu  1996). 

In  the  active  system,  an  integrated  thermal  signal  to 
the  thermoregulatory  controller  is  composed  of  the 
weighted  thermal  input  from  thermal  receptors  at  vari¬ 
ous  sites  distributed  throughout  the  body.  The  inte¬ 
grated  body  temperature  is  weighted  using  the  core, 
muscle  and  skin  compartment  temperatures.  The  affer¬ 
ent  signal  is  the  difference  between  this  temperature 
and  its  threshold,  which  activates  thermoregulatory 
mechanisms  including  vasomotor  changes,  sweat  pro¬ 
duction  and  metabolic  heat  production  (Xu  and  Wer¬ 
ner  1997).  Shivering  thermogenesis  (i.e.  part  of 
metabolic  heat  production)  is  a  function  of  core  and 
mean  skin  temperatures,  and  includes  an  intensity 
adjustment,  maximal  capabUity,  shivering  exhaustion, 
and  inhibition  due  to  a  low  core  temperature  (Xu  et  al. 
2005).  The  maximal  shivering  intensity  was  estimated 
from  the  height,  weight,  V02^ax  ^^d  age  (Eyolfson 
et  al.  2001). 

CTM  inputs  include  individual  characteristics  (i.e. 
height,  weight,  percentage  of  body  fat,  age, 
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and  exercise  intensity,  as  well  as  environmental 
(i.e.  temperature,  humidity,  and  wind  velocity)  and 
clothing  (clothing  insulation  do,  moisture  permeability 
index  („,)  properties  for  each  of  the  six  cylinders. 

Experimental  data 

Ten  male  subjects  volunteered  for  the  experiment. 
Their  physical  characteristics  are  listed  in  Table  1.  This 
study  was  approved  by  the  appropriate  Scientific  and 
Human  Use  Review  Boards  of  the  U.S.  Army 
Research  Institute  of  Environmental  Medicine  and  US 
Army  Medical  Research  and  Materiel  Command.  The 
subjects  volunteered  after  being  fully  informed  of  the 
requirements  and  risks  associated  with  the  research. 
For  analysis  purpose,  subjects  were  grouped  according 
to  their  percentage  of  body  fat  as  low  body  fat  (LF 
<14.7%,  five  subjects  A— E)  or  normal  body  fat  (NF 
>15.0%,  five  subjects  F-J).  The  percent  body  fat  was 
determined  from  the  equation  of  Durnin  and  Womers- 
ley  using  four  skinfolds  (Durnin  and  Womersley  1974). 
Subjects  walked  at  0.44  m  s~^,  at  two  different  immer¬ 
sion  depths  (waist  to  the  iliac  crest  and  chest  to  the  nip¬ 
ple)  and  at  two  different  water  temperatures  (10  and 
15°C).  The  air  temperature  near  the  head  and  chest 
ranged  from  15  to  22°C,  depending  on  the  water  tem¬ 
perature.  Subjects  put  their  arms  on  a  plastic  platform 
to  keep  their  arms  out  of  water.  AU  trials  were  sepa¬ 
rated  by  at  least  one  day  of  rest.  Subjects  wore  the  US 
Army  battle-dress  uniform  (BDU;  dry  clothing  insula¬ 
tion  value  of  1.3  do),  socks,  and  neoprene  water  shoes 
during  exercise-cold  water  immersion.  No  caffeine  or 
alcohol  was  consumed  on  the  test  day.  Rectal  tempera¬ 
ture  (r<, )  was  measured  using  a  rectal  probe  (YSI,  Yel¬ 
low  Springs,  OH,  USA)  inserted  10  cm  beyond  the 
anal  sphincter.  Skin  temperature  (°C)  and  heat  flow 
(W  m“^)  were  measured  by  heat  flow  sensors  with  an 
integrated  thermistor  (Concept  Engineering,  Old  Say- 
brook,  CT,  USA)  attached  on  the  skin  surface  at  eight 
sites  (on  the  right  side):  anterior  aspect  of  forearm, 
forehead,  subscapular,  triceps,  pectorahs  major,  abdo¬ 
men  (7.5  cm  lateral  to  umbihcus),  anterior  thigh,  and 
calf.  Temperature  data  were  collected  every  15  s  dur¬ 
ing  treadmill  walking  (PX  1006,  National  Instruments). 
Oxygen  uptake  (UO2)  was  determined  using  an  online 
metabolic  analysis  system  (Sensormedics  Vinax,  Yorba 
Linda,  CA)  before  exercise  and  ~every  20  min  during 
walking.  Subject’s  expired  air  was  collected  for  5  min 
each  time  using  a  mouthpiece  and  noseclip.  The  order 
for  the  four  trials/test  conditions  was  randomized  and 
counterbalanced.  Each  subject’s  trial  began  at  approxi¬ 
mately  the  same  time  each  day.  They  walked  on  an 
underwater  treadmill  until  one  of  the  following 


occurred:  their  rectal  temperature  reached  SS-S'C,  they 
had  exercised  for  4  h,  the  subject  asked  to  stop,  or  sub¬ 
ject’s  exposure  was  stopped  by  the  principal  investiga¬ 
tor.  Each  combination  was  to  be  tested  once,  with 
seven  of  ten  subjects  completing  four  trials.  Two  sub¬ 
jects  completed  three  trials  and  one  subject  completed 
only  two  trials,  they  stopped  their  participation  and 
voluntarily  withdrew  from  the  study.  Thus,  the  total 
number  of  trials  completed  was  36. 

CTM  input  adjustment  for  partial  immersion 

From  a  simulation  viewpoint,  the  environments  that 
subjects  were  exposed  to  were  complex  as  conditions 
were  not  uniform,  i.e.  subjects  were  exposed  to  both  air 
and  water  at  the  same  time.  The  clothing  was  either 
wet  or  becoming  wet  from  the  absorption  of  water  dur¬ 
ing  the  experiments.  Of  the  six  CTM  cylinders,  the 
head,  arms  and  hands  were  always  exposed  to  the  air, 
and  the  legs  and  feet  were  always  immersed.  The  trunk 
cylinder  was  exposed  to  both  air  and  water  of  varying 
depth  and  the  clothing  insulation  covering  the  trunk 
was  in  transition.  Thus,  the  environmental  and  clothing 
parameters  for  the  trunk  needed  to  be  modified  to  take 
these  factors  into  account.  Based  on  the  experimental 
design,  it  was  assumed  that  20%  of  the  torso  was 
exposed  to  air  while  80%  was  immersed  during  CH, 
whereas  during  WA,  80%  of  the  torso  was  exposed  to 
air  while  20%  was  immersed.  Therefore,  the  adjusted 
parameters  were  calculated  using  area  weighted  values 
for  air  and  water  which  reflected  the  immersion  depth. 

Simulation  and  evaluation 

CTM  simulations  were  run  for  36  individual  trials  in 
CHIO,  CH15,  WAIO  and  WA15.  The  measured  was 
compared  with  the  predicted  T^,  using  paired  t  tests 
and  the  root  mean  square  deviation  (RMSD).  RMSD 
quantifies  the  average  difference  between  predicted 


Table  1  Individual  physical  characteristic  of  subjects 


Subject 

Height 

(m) 

Weight 

(kg) 

Body 

fat(%) 

Age 

F02max 
(ml  kg“'  min“^) 

A 

1.78 

86.6 

14.7 

19 

48.2 

B 

1.93 

77.3 

8.1 

25 

38.9 

C 

1.83 

70.7 

14 

20 

44.4 

D 

1.80 

72.4 

14.7 

19 

50-4 

E 

1.73 

64.1 

12.9 

18 

47.2 

F 

1.75 

74.5 

17 

20 

42.5 

G 

1.73 

66.3 

15 

18 

47.8 

H 

1.70 

65.4 

19 

20 

47.2 

I 

1.75 

76.1 

16.4 

18 

48.5 

J 

1,77 

80.4 

17.7 

19 

57.1 
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and  observed  measurement  across  time  (Haslam  and 
Parsons  1994;  Kraning  and  Gonzalez  1997).  The 
RMSD  was  calculated  between  the  observed  and  pre¬ 
dicted  temperatures  using  a  10-min  time  interval.  Data 
were  analyzed  using  a  two-factor  (trial  x  time) 
repeated-measures  ANOVA  to  compare  between  pre¬ 
dicted  and  observed  values.  Trials  are  defined  as  the 
measured  and  predicted  values.  When  significant  F 
ratios  were  calculated,  paired  comparisons  were  made 
post  hoc  using  a  Tukey  HSD  test.  The  level  of  signifi¬ 
cance  for  differences  reported  is  P  <  0.05.  To  further 
investigate  the  potential  influences  of  shivering  ther¬ 
mogenesis  and  muscle  blood  flow  in  legs  on 
responses  in  the  LF  group,  CTM  was  modified  slightly 
and  run  under  two  particular  conditions.  The  modifica¬ 
tions  included  inhibition  of  shivering  or  an  increase  in 
basal  muscle  blood  flow  rate  in  legs  from  1.6  ml/min 
per  100  ml  tissue  (Xu  and  Werner  1997)  to  5  ml/min 
per  100  ml  tissue  (Johnson  1989).  Limbs  and  extremi¬ 
ties  have  large  vasomotor  regulation  (Toner  and 
McArdle  1996)  and  thus,  blood  flow  in  these  areas  may 
be  high  either  before  or  during  exercise  in  cold  water 
and  therefore  play  a  critical  role  in  T^. 

Results 

Figure  1  shows  the  measured  of  each  individual  in 
NF  group  and  mean  of  the  predicted  at  CHIO, 
CH15,  WAIO  and  WA15.  The  observed  increased 
slightly  during  the  beginning  of  immersion  due  to 
vasoconstriction,  and  then  started  to  fall  after  reaching 
its  peak  in  about  20  min.  The  predicted  was  close  to 
the  measured  7).  and  within  the  span  of  experimental 
values.  Paired  t  tests  indicated  that  differences 
between  the  measured  and  predicted  were  not  sig¬ 
nificant  (P  >  0.05)  at  CHIO,  CH15  and  WAIO,  but  at 
WA15  there  were  no  significant  differences  except  for 
the  data  at  110  min  (P<  0.0001  or  0.0001).  Figure  2 
shows  the  measured  of  each  individual  in  LF  group 
and  mean  of  the  predicted  for  each  of  the  trials. 
responses  of  the  LF  group  appeared  to  be  different 
from  the  NF  group.  At  the  beginning  of  immersion, 
the  Pg  of  LF  group  increased  for  a  shorter  period,  and 
then  started  to  drop  quickly.  For  one  subject, 
dropped  immediately  after  immersion.  Paired  t  tests 
indicated  that  the  differences  between  measured 
and  predicted  were  significant  (P  <  0.0001)  at  all 
four  conditions.  Figure  2  also  shows  the  mean  of  the 
adjusted  predicted  T^,  i.e.  the  prediction  with  an 
increased  muscle  blood  flow  in  legs;  the  adjusted  pre¬ 
diction  results  match  more  closely  with  observed  val¬ 
ues  than  the  original  prediction. 


Table  2  shows  the  mean  (SD)  of  the  measured  and 
predicted  at  the  end  of  each  trial,  the  RMSD  and  the 
average  (SD)  immersion  time  for  each  trial.  The  pre¬ 
dicted  P(,  for  the  NF  group  was  within  the  range  of  the 
measured  ±  SD  in  all  four  trials,  with  RMSD  values 
ranging  from  0.15  to  0.32“C.  These  RMSD  values  were 
within  or  close  to  the  SD  values  of  NF.  For  the  LF 
group,  the  predicted  was  vrithin  the  range  of  the 
measured  ±  SD  only  at  CH15  and  was  higher  than 
the  measured  for  the  other  three  trials  with  the 
RMSD  values  ranging  from  0.3  to  0.6°C.  When  the 
muscle  blood  flow  was  adjusted  for  LF  group,  the  pre¬ 
dicted  Pp  at  the  end  of  immersion  were  35.43°C  for 
CHIO,  35.13°C  for  CH15,  35.67°C  for  WAIO  and 
35.65°C  for  W15  separately  and  their  correspondent 
RMSD  were  0.28,  0.26,  0.24  and  0.26°C,  respectively. 
The  adjusted  predicted  P^  was  closer  to  observed 
and  RMSD  were  smaller  in  comparison  with  the  origi¬ 
nal  RMSD.  The  immersion  time  for  the  NF  group  was 
about  75-123%  longer  than  the  LF  group. 

The  effects  of  shivering  intensity  on  responses  are 
shown  in  Fig.  3,  the  results  for  Subject  C  at  CHIO  and 
Subject  E  at  WAIO  are  depicted  as  an  example.  When 
shivering  was  inhibited  in  the  CTM  simulation,  the  pre¬ 
dicted  Pj.  was  lower  than  the  original  simulation  when 
shivering  was  not  inhibited.  Thus,  the  predicted  P^  was 
closer  to  the  measured  P^  for  both  subjects  at  CHIO  or 
WAIO  conditions  and  the  RMSD  was  reduced  from 
0.38  to  0.24'’C  for  Subject  C  and  from  0.59  to  0.47°C  for 
Subject  E  when  shivering  was  inhibited  in  the  simula¬ 
tion.  When  the  muscle  blood  flow  in  legs  was  increased 
to  5.0  mFmin  per  100  ml  tissue,  the  predicted  P^  was 
improved  significantly  and  matched  the  measured 
well.  After  adjusting  for  muscle  blood  flow,  the 
RMSDs  were  0.14°C  for  Subject  C  at  CHIO  and  0.15 
for  Subject  E  at  WAIO. 


Discussion 

The  CTM  was  successfully  used  to  simulate  individual 
thermal  responses  to  cold  exposure.  Only  a  few  models 
consider  individual  characteristics  and  aU  were 
designed  for  hot  or  warm/cool  environments.  Havenith 
(2001)  developed  an  individualized  model  of  human 
thermoregulation  for  simulation  of  heat  stress 
responses  while  Zhang  and  her  colleagues  (2001) 
developed  a  body  builder  model  and  incorporated  this 
into  the  UC  Berkeley  thermal  comfort  model  to  pre¬ 
dict  variances  in  thermal  response  between  individuals. 
The  CTM  was  designed  for  the  cold  with  the  capability 
for  individualization  and  provides  a  significant  upgrade 
to  the  original  model  and  could  be  utilized  for  many 
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Fig.  1  Time  versus  individual  measured  core  temperatures  (each  15°C  (CH15),  Waist  10°C  (WAIO)  and  Waist  IS^C  (WA15). 
symbol  represents  one  individual)  and  mean  of  predicted  core  Asterisk  denotes  predicted  value  is  significantly  different  than  the 
temperature  {line)  in  subjects  with  a  %  fat  of  15— 19%  (NF).  Trials  mean  measured  value  at  specified  time  points 

consisted  of  walking  at  0.44  m  s“^  at  chest  10°C  (CHIO),  Chest 


WA10NF  WA15NF 


applications.  Significant  inter-individual  variance  in 
rectal  temperature,  skin  temperature  and  energy 
expenditure  in  response  to  mild  cold  has  been  reported 
(Marken  Lichtenbelt  et  al.  2002)  and  data  from  the 
current  study  clearly  demonstrate  inter-individual 
differences  during  light  exercise  in  cold  water  (e.g.,  one 
LF  subject  reached  of  35.5°C  in  about  20  min  at 
CHIO  conditions  while  two  NF  subjects  maintained 
their  above  35.5°C  after  120  min).  This  study  found 
that,  in  these  ten  subjects  during  immersion  in  cold 
water  at  two  depths,  CTM  predicted  responses  with 
acceptable  accuracy  for  certain  individuals,  mostly  NF 
people,  and  tended  to  overpredict  responses  for 
other  individuals,  mostly  LF  people. 

CTM  requires  height,  weight,  percentage  of  body 
fat,  age,  and  maximal  oxygen  consumption  to  predict 


individual  thermal  responses.  Height,  weight  and  per¬ 
centage  of  body  fat  determine  the  geometrical  sizes  of 
each  cylinder  and  layer  (Xu  1996).  Body  fat,  age  and 
maximal  oxygen  consumption  are  used  to  determine 
the  shivering  intensity,  maximal  shivering  capacity  and 
shivering  fatigue  (Xu  et  al.  2005).  In  a  previous  study 
(Xu  et  al.  2005),  CTM  predictions  were  in  good  agree¬ 
ment  with  measured  observations  for  ten  male  and 
female  subjects  immersed  in  8-10“C  water  who  had 
body  fat  percentages  ranging  from  14.9  to  33.6%.  The 
current  study  further  demonstrated  that  CTM  predic¬ 
tion  for  the  NF  group  was  in  good  agreement  with 
observation  during  partial  immersion  at  two  depths 
during  slow  walking.  The  five  inputs  appear  to  ade¬ 
quately  describe  individual  responses  to  cold-stress  for 
a  certain  group  of  individuals,  i.e.  people  with  a  normal 
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Fig.  2  Time  versus  individual  measured  core  temperatures  (each 
symbol  represents  one  individual),  mean  of  predicted  core  tem¬ 
perature  (line)  in  subjects  with  a  %  fat  of  <14.7%  (LF)  and  mean 
of  predicted  core  temperature  with  increased  muscle  blood  flow 
in  legs  (adjusted  predicted  core  temperature,  thin  line).  Trials 


consisted  of  walking  at  0.44  m  s“^  at  chest  10°C  (CHIO),  chest 
15°C  (CH15),  waist  10°C  (WAIO)  and  waist  15°C  (WA15).  Aster¬ 
isk  denotes  predicted  value  is  significantly  different  than  the  mean 
measured  value  at  specified  time  points 


percentage  of  body  fat.  What  other  physiological  infor¬ 
mation,  inputs  or  anthropometric  data  (e.g.  length  of 
the  leg)  would  be  required  to  describe  individual 
responses  to  cold,  especially  in  individuals  with  a  low 
percentage  of  body  fat? 

There  are  initial  differences  before  cold  exposure 
between  the  measured  and  predicted  T^.  Predicted  7). 
was  fixed  at  a  constant  value  of  36.8°C  in  the  current 
CTM  construct,  whUe  the  measured  pre-exposure  7), 
varied  between  and  within  subjects.  As  shown  in  Fig.  2, 
the  initial  temperatures  of  the  LF  group  ranged  from 
36.2  to  36.9'’C  in  trial  W15.  Subject  C  (symbol  x)  had 
an  initial  temperature  of  36.2°C  in  W15  and  36.7°C  in 
CH15.  As  the  endurance  times  for  the  LF  group  were 
short  due  to  rapid  core  cooling,  initial  increased  in 


relative  importance.  With  a  given  cooling  rate,  a  lower 
initial  caused  subjects  to  reach  the  experimental 
end-point  in  a  shorter  time.  Figure  2  shows  that  the 
measured  for  most  subjects  was  lower  than  the  pre¬ 
dicted  T^.  Individuals  with  low  initial  core  tempera¬ 
tures  and  low  percent  body  fat  appear  more  vulnerable 
to  h)q3othermia.  It  might  improve  prediction  accuracy 
if  CTM  was  modified  to  allow  the  input  of  an  individ¬ 
ual’s  initial  7).. 

CTM  considered  two  basic  roles  of  body  fat  in  main¬ 
taining  body  heat  balance  and  core  temperature  during 
cold  exposure:  insulation  and  attenuation  of  shivering 
heat  production.  Fat  provides  insulation,  due  to  the 
low  heat  conductivity,  reported  to  be  0.15-0.5  W/ 
m°C  (Werner  and  Buse  1988),  thus  providing  higher 
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Table  2  Mean  (SD)  measured  and  predicted  core  temperature 
(rj  at  the  end  of  each  trial,  root  mean  square  deviation  (RMSD), 
and  mean  (SD)  immersion  time  during  each  trial  for  the  normal 
fat  (NF)  and  low  fat  (LF)  groups 


Trials 

T;  (°C)  measured 
predicted 

RMSD 

(“C) 

Time 

(min) 

NF 

CHIO 

35.85  (0.22) 

35.72 

0.32 

96(35) 

CH15 

36.08  (0.35) 

35.85 

0.23 

140  (19) 

WAIO 

36.23  (0.25) 

35.98 

0.15 

123  (10) 

WA15 

36.04  (0.29) 

36.26 

0.30 

172  (12) 

LF 

CHIO 

35.54  (0.05) 

36.13 

0.52^ 

43  (19) 

CH15 

35.87  (0.37) 

36.14 

0.31® 

80  (32) 

WAIO 

35.60  (0.18) 

36.56 

0.60® 

54  (23) 

WA15 

35.65  (0.07) 

36.46 

0.58® 

98  (57) 

“  RMSD  with  adjusted  prediction  are  0.28, 0.26, 0,25,  and  0.26°C, 
respectively 


thermal  resistances  than  either  skin  or  muscle  tissue 
(Toner  and  McArdle  1988).  The  mean  thickness  of 
subcutaneous  fat  in  male  subjects  with  8.8  to  27.2% 
percent  body  fat  ranged  from  1.1  mm  to  9.6  mm  (Veic- 
steinas  et  al.  1982).  That,  with  a  given  temperature 
difference  over  the  fat  layer  and  a  constant  heat  con¬ 
ductivity,  would  represent  about  an  8.7-fold  difference 
in  heat  flow,  i.e.  heat  flow  with  8.8%  body  fat  is  about 
8.7  times  as  large  as  the  heat  flow  with  27.2%  body  fat. 
Based  on  the  current  CTM  construct,  the  fat  thickness 
is  ~3.7  mm  for  10%  body  fat  and  9.6  mm  for  25%  in  a 
human  with  height  of  1.83  m  and  weight  of  70.7  kg. 
Many  studies  have  clearly  demonstrated  the  benefits  of 
body  fat  in  maintenance  of  core  temperature  during 
cold-water  immersion.  A  relatively  fat  person  with 
large  subcutaneous  fat  and  minimal  perfusion  of  the 
skin  and  extremities  can  withstand  the  challenges  of 
cold  water  more  effectively  than  a  leaner  counterpart 
(Toner  and  McArdle  1996;  Stocks  et  al.  2004).  The  role 
of  fat  as  insulation  becomes  even  more  important  dur¬ 
ing  exercise,  as  the  insulation  provided  by  muscle, 
which  provides  60-75%  of  body  insulation,  is 
decreased  (Veicsteinas  etal.  1982;  Park  etal.  1984). 
However,  if  shivering  can  increase  heat  production 
substantially,  the  insulation  provided  by  fat  may  be  less 
important.  During  90  min  immersion  in  18,  22  and 
26°C  water,  low  fat  male  subjects  were  able  to  maintain 
a  similar  when  compared  to  their  high  fat  counter¬ 
parts  due  to  greater  shivering  thermogenesis  (GHck- 
man-Weiss  et  al.  1991).  In  addition,  shivering 
metabolism  is  inversely  proportional  to  the  percentage 
of  body  fat  during  cold-water  immersion  (Tikuisis  et  al. 
1988).  According  to  the  CTM  construct,  low  fat  indi¬ 
viduals  have  a  thin  fat  layer  and  a  high  metabolic  rate. 
Yet  the  high-predicted  metabolic  rate  does  not  appear 


Tc{°C)  Subject  C  cm  0 


Time  (min) 


Fig.  3  Measured  and  predicted  core  temperatures  for  subject  C 
at  chest  10‘’C  (CHIO)  immersion  {top)  and  for  subject  E  at  waist 
10°C  (WAIO)  immersion  {bottom),  predicted  core  temperatures 
while  shivering  was  inhibited  or  leg  muscle  blood  flow  was  in¬ 
creased  to  5.0  ml/min  per  100  ml  tissue  in  the  CTM  simulation 

to  be  the  reason  that  the  predicted  7),  is  higher  than  the 
measured  T^.  As  shown  by  the  two  individual  examples 
in  Fig.  3,  even  when  shivering  heat  production  is  inhib¬ 
ited  (i.e.  set  to  zero)  in  CTM,  the  predicted  is  still 
higher  than  the  measured  for  both  subject  C  in  the 
CHIO  trial  and  subject  E  in  the  WAIO  trial.  Therefore, 
the  insulation  function  of  body  fat  and  shivering  heat 
production  appear  not  to  explain  the  phenomena,  sug¬ 
gesting  other  mechanisms  may  contribute  to  the  rapid 
decline  in  observed  m  the  LF  group. 

Blood  flow  influences  thermal  responses  to  cold  but 
its  effects  have  not  been  fuUy  characterized  yet,  as  the 
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measurement  of  blood  flow,  especially  muscle  blood 
flow,  is  difficult.  During  cold  water  immersion,  the  high 
gradient  between  the  core  and  muscle  in  the  extremi¬ 
ties  [e.g.  18°C  during  8°C  water  immersion  (Bristow 
etal.  1994)]  and  high  amount  of  muscle  mass,  i.e. 
~40%  of  body  mass  could  facilitate  heat  exchange 
between  the  core  and  extremity  even  with  relatively 
low  blood  flows  in  muscle.  The  limbs  and  extremities, 
body  areas  that  have  large  vasomotor  regulation,  have 
two  adjustable  thermal  resistances:  muscle  blood  flow 
and  skin  blood  flow  (Toner  and  McArdle  1996).  At  the 
beginning  of  cold  immersion,  increased  slightly  for 
~10  min  with  LF  group  and  ~20  min  with  NF  group 
due  to  vasoconstriction  and  then  started  to  decrease. 
As  shown  in  Fig.  2,  when  the  leg  muscle  blood  flow  (i.e. 
blood  flow  before  cold  exposme)  is  increased  to  5  ml/ 
min  per  100  ml  tissue  (from  1.6  ml/min  per  100  ml  tis¬ 
sue),  the  predicted  for  LF  agreed  with  the  measured 
more  closely  in  aU  four  trials.  Figure  3  indicated  that 
is  more  sensitive  to  changes  in  leg  muscle  blood  flow 
than  to  the  changes  in  shivering.  Hence,  the  CTM  sim¬ 
ulations  support  an  interpretation  that  high  leg  muscle 
blood  flow  causes  a  rapid  decline  in  T^.  A  recent  study 
demonstrated  that  individual  variance  in  muscle  blood 
flow  is  greater  than  individual  variance  in  skin  blood 
flow  (Kenny  et  al.  2006).  Thus  it  is  reasonable  to  postu¬ 
late  that  individual  differences  in  muscle  blood  flow 
partly  causes  the  individual  thermal  responses  to  par¬ 
tial  immersion.  Furthermore,  body  fat  could  modulate 
the  effects  of  muscle  blood  flow  on  response.  First,  a 
low  percent  body  fat  is  associated  with  a  thin  fat  layer 
and  the  muscle  is  close  to  the  skin  surface,  which  could 
cause  the  outer  layer  of  muscle  in  legs  to  cool  down 
quickly  at  the  beginning  of  immersion,  as  well  as  cool 
down  the  temperature  of  venous  blood,  leading  to  core 
coohng.  This  seems  consistent  with  a  previous  observa¬ 
tion  that  older  swimmers  and  boys  of  given  trunk  fat 
thickness  generally  had  less  hmb  fat  and  cooled  faster 
than  younger  swimmers  and  girls  (Sloan  and  Keatinge 
1973).  Second,  body  fat  composition  may  also  affect 
blood  flow  rates  and  vasomotor  regulations.  Lean  indi¬ 
viduals  (ectomorphs)  have  a  higher  blood  volume  per 
unit  weight  than  fat  individuals  (endomorphs)  (Gre- 
gersen  and  Nickerson  1950).  Blood  flow  in  muscle  is 
more  severely  restricted  in  fat  subjects  than  m  thin  sub¬ 
jects  in  water  of  critical  temperature  (Park  et  al.  1984). 
This  might  partly  explain  the  differences  in  thermal 
responses  of  the  NF  and  LF  groups  to  cold  immersion. 

CTM  simulated  effects  of  the  muscle  blood  flow  in 
legs  on  Tj.  cooling  and  the  prediction  indicated  that  leg 
muscle  blood  flow  in  legs  is  a  possible  contributor  to 
the  rapid  decline  in  7),  for  LF  people.  In  the  experi¬ 
ments,  the  subjects  were  walking  at  a  speed  of  0.44  m/s 


on  a  treadmill  during  immersion  and  exercise  increases 
muscle  blood  flow  rates  in  legs.  Other  reports  indicate 
blood  flow  in  the  calf  during  rest  of  ~3.6  ml/min  100  ml 
tissue  (Savard  et  al.  1985)  and  blood  flow  to  inactive 
muscle  of  ~5-10  ml/min  per  100  g  muscle  (Johnson 
1989).  Thus  the  revised  muscle  blood  flow  rates  in  the 
CTM  simulation,  5.01/min  per  100  ml  tissue,  seem 
physiologically  reasonable  and  the  adjusted  predicted 
consistent  with  observations.  Exercise  substantially 
lowers  body  insulation  during  cold  immersion,  proba¬ 
bly  by  increasing  muscle  blood  flow  (Sloan  and  Keat¬ 
inge  1973;  Yeon  et  al.  1987).  The  contribution  from 
cooled  limb  blood  flow  to  a  decrease  in  has  also 
been  observed  during  cold-water  immersion  where 
blood  flow  to  the  Umbs  was  occluded  (Mittleman  and 
Mekjavic  1988).  Even  under  neutral  environmental 
conditions,  convective  heat  transfer  through  blood  flow 
likely  contributes  to  core  temperature  changes  and 
heat  redistribution  between  the  muscle  and  core  region 
during  the  postexercise  period  (Kenny  etal.  2006). 
Due  to  the  possible  contribution  of  leg  muscle  blood 
flow  to  the  observed  responses,  individuals  with  low 
body  fat  percentages  are  potentially  at  greater  risk  dur¬ 
ing  exercise  while  only  partially  immersed. 

One  limitation  of  the  present  study  was  the  separa¬ 
tion  of  the  groups  for  analysis.  The  LF  group  had  only 
one  very  low  fat  subject  (i.e.  8.1%)  and  the  NF  group 
had  one  high  fat  subject  (i.e.19%),  but  the  fat  differ¬ 
ences  between  the  two  groups  were  not  large.  There¬ 
fore,  further  study  is  needed  to  examine  whether  the 
relation  between  CTM  accuracy  and  fat  %  observed  in 
this  study  could  be  expanded  to  a  broader  scope.  How¬ 
ever,  findings  of  this  study  do  suggest  that  biophysical 
and  physiological  parameters  (i.e.  height,  weigh,  and 
body  fat  %  and  shivering  heat  production)  do  not  fully 
explain  large  individual  differences  in  7).  responses  to 
partial  immersion.  Muscle  blood  flow,  either  basal  val¬ 
ues  or  increases  due  to  exercise,  and  blood  flow  redis¬ 
tribution  could  play  critical  roles  in  thermal 
equilibrium  and  account  for  individual  differences  in 
responses  during  immersion.  Two  prior  studies  demon¬ 
strate  the  need  for  an  independent  compartment,  i.e. 
muscle,  in  addition  to  core  and  shell,  to  assess  the  heat 
balance  status  of  the  body  (Tflcuisis  2003;  Jay  etal. 
2006).  Convection  through  blood  flow  would  be  an 
effective  way  to  exchange  heat  among  these  three  com¬ 
partments.  Although  thermal  responses  to  cold  have 
been  well  studied,  there  is  limited  quantitative  infor¬ 
mation  and  knowledge  about  how  cold  impacts 
regional  blood  flows.  Modeling  techniques  have  been 
long  used  to  investigate  the  effects  of  blood  flow  on  local 
tissue  temperature  distribution  during  induced  regional 
hyperthermia  therapies  (Brinck  and  Werner  1995; 
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Kolios  et  al.  1998;  Alekssev  et  al.  2005).  CTM  simula¬ 
tion  demonstrated  the  potential  contribution  of  muscle 
blood  flow  to  cooling  by  means  of  a  whole  body 
thermoregulatory  model. 

Immersion  at  various  depths  represents  a  relatively 
complex  set  of  environmental  conditions.  Heat  loss  in 
air  is  significantly  different  from  heat  losses  in  water. 
Clothing  absorbs  water  through  capillary  action  and 
becomes  progressively  wetter,  thus  the  thermal  insula¬ 
tion  values  change  over  time.  According  to  the  basic 
CTM  construct,  each  cylinder  should  be  exposed  to  a 
uniform  environment  with  a  specific  set  of  physical  and 
clothing  parameters.  To  adapt  the  model  for  partial 
immersion,  adjustments  based  on  the  percentage  of  the 
cylinder  that  was  immersed  were  made  to  the  inputs 
for  the  torso  section.  If  the  percentage  was  changed  by 
±5.0%  from  the  80.0%  during  CH,  the  deviations  of 
CTM  prediction  at  end  of  180  min  immersion  were  less 
than  Q.1°C.  The  CTM  predictions  neglected  the  influ¬ 
ence  of  exercise  for  simplicity,  as  the  exercise  load  (i.e. 
walking  at  0.44  m/s)  was  low.  In  addition,  for  modeling 
the  trials,  measured  heat  transfer  coefficients  during 
immersion  were  used,  instead  of  model  default  values, 
to  make  the  simulation  more  reahstic  and  account  for 
the  effect  of  water  movement  on  convective  heat  loss. 
Results  from  this  study  indicate  that  adjustments  (i.e. 
parameters  for  the  torso  section)  to  the  model  based 
on  best  available  information  are  necessary  to  ensure 
that  the  inputs  and  predictions  accurately  represent  the 
actual  conditions,  especially  when  CTM  is  used  to  pre¬ 
dict  thermal  responses  in  the  field.  To  ensure  that  the 
inputs  are  realistic,  rigorous  evaluation  to  identify  aU 
possible  scenarios  is  required.  When  the  information 
available  is  limited,  it  might  be  necessary  to  run  ther¬ 
mal  response  simulations  for  the  worst  and  best  cases. 
The  results,  which  indicate  that  the  basic  CTM  can  be 
successfully  adapted  to  the  more  complex  scenario  of 
partial  immersion,  suggest  that  CTM  could  be  applied 
to  other  partial  immersion  scenarios  such  as  one  foot 
or  one  leg  immersion. 


Conclusions 

CTM  predicted,  in  this  study  with  acceptable  accuracy, 
core  temperature  responses  to  immersion  at  various 
depths  for  certain  individuals,  i.e.  those  with  normal 
body  fat  compositions.  The  differences  in  the  measured 
and  predicted  core  temperatures  for  people  with  low 
body  fat  composition  are  Ukely  due  to  initial  tempera¬ 
ture  (pre-exposure  temperature)  and  muscle  blood 
flow  in  the  legs.  Simulation  analysis  suggested  that  the 
leg  muscle  blood  flow  may  contribute  significantly  to 


the  rapid  drop  in  core  temperature  and  partly  account 
for  individual  thermal  responses  to  partial  immersion. 
CTM  can  be  adapted  to  non-uniform  environments  by 
adjusting  the  input  accordingly.  When  scenarios  are 
complex  and/or  available  information  is  minimal,  it 
would  be  necessary  to  run  CTM  for  the  worst  and  best 
cases,  and  the  influential  factors,  e.g.  initial  tempera¬ 
ture  and  leg  muscle  blood  flow,  should  be  considered. 
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